Higher Order Slow-Roll Predictions for Inflation 
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We study the WKB approximation beyond leading order for cosmological perturbations during 
inflation. To first order in the slow-roll parameters, we show that an improved WKB approximation 
leads to analytical results agreeing to within 0.1% with the standard slow-roll results. Moreover, 
the leading WKB approximation to second order in the slow-roll parameters leads to analytical 
predictions in qualitative agreement with those obtained by the Green's function method. 
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PACS numbers: 98.80.Cq, 98.80.-k 

Introduction. It is nowadays common to state that we 
are in the era of precision cosmology. The implication of 
present and future data sets will be able to discriminate 
among different inflationary models 0. For this reason, 
the comparison of inflationary models with observations 
requires theoretical advances in the predictions for the 
power spectrum of primordial perturbations beyond the 
lowest order in the slow-roll parameters worked out by 
Stewart and Lyth 0. Such slow-roll parameters quan- 
tify the deviation from an exactly exponential expansion 
during inflation and are related for a canonical inflaton to 
the derivative of the potential. Among the many different 
conventions, we find it convenient to use the hierarchy of 
horizon flow functions e, 0, defined as e i+ i = £{/ (Ha), 
with ei = —H/H 2 , H = a/ a the Hubble parameter, a 
the scale factor and dots denote derivatives respect to the 
cosmic time. 

The search for deviations from a simple power-law pa- 
rameterization of the primordial power spectrum in the 
first year of WMAP data has begun Q, with however 
no statistically conclusive evidence for a significant de- 
viation. Stronger conclusions about any deviations will 
be made possible by the better resolution of the Planck 
satellite J5|- In the context of inflation, the deviations 
from power-law spectra considered so far correspond to 
some of the predictions beyond the first order in the slow- 
roll parameters and are expected in general. 

The purpose of this letter is to show how the primordial 
power spectrum of cosmological perturbations generated 
in a single field inflationary model can be analytically pre- 
dicted by the WKB method with sufficient accuracy. We 
shall also show how two different analytical approxima- 
tions, such as the Green's function method (GFM hence- 
forth) [3 and the WKB method 0,|| used here, agree to 
second order on the polynomial structure of the results 
in the horizon flow functions during inflation, but differ 
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in the numerical coefficients of the second order terms 
(the same polynomial structure in the s pect ral indices is 
found by the uniform approximation yj, UJJ, III! ) • 
Cosmological perturbations. Let us begin by recalling 
that scalar (density) and tensor (gravitational wave) fluc- 
tuations on a Robertson- Walker background are given 
respectively by \i = /is = a Q (Q is the Mukhanov vari- 
able and fi = fiT = ah (his the amplitude of the two 
polarizations of gravitational waves 

[D|[lil). The func- 
tions n must satisfy the one-dimensional Schrodinger-like 
equation 



H = 



together with the initial condition 



lim /j,(k,rf) 



V2k 



(1) 



(2) 



In the above rj is the conformal time, fc is the wave- 
number, and 



ft 2 (fc,77) = fc 2 - — 



(3) 



where 



2 4>' I H for scalar and z — zj- = a for 
tensor perturbations (4> is the homogenous inflaton and 
primes denote derivatives with respect to rj) . The dimen- 
sionless power spectra of scalar and tensor fluctuations 
are then given by 



fc 3 



4fc 3 



fir 



(4a) 



and the spectral indices and runnings by 
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dlnP c 
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k=fc. 
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k=k. 



where fc* is an arbitrary pivot scale. We also define the 
tcnsor-to-scalar ratio 



R = 



Vh 



(4d) 



k=k. 
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Slow-roll and WKB approximations. One of the very few 
cases for which the equations for cosmological pertur- 
bations can be integrated exactly is that of power-law 
inflation [TH ITt^ , where the inflaton is a canonical scalar 
field with an exponential potential. In such a case, t\ 
is constant in time and the hierarchy of the horizon flow 
functions is therefore truncated with ti = for i > 2. The 
original slow-roll approximation corresponds to consider- 
ing both ei and £2 constant in time for any potential 0. 

Although apparently inappropriate for studying phys- 
ical problems such as the hydrogen atom [l7| (or cos- 
mological perturbations 0), the WKB method can be 
cleverly applied after suitable redefinitions of the wave- 
function (for Fourier modes k) and variables (with the 
corresponding "frequency" Q replaced by a new expres- 
sion uj, as given in detail in Refs. 0,0)- This improved 
WKB method applied to cosmological perturbations with 
a linear turning point in to does not however predict am- 
plitudes with a sufficient accuracy to lowest order . We 
have shown in Ref. how the prediction for the ampli- 
tude may be improved by using a next-to-leading WKB 
approximation. Our method involved an adiabatic ex- 
pansion and the numerical evaluation of the higher-order 
coefficients which yielded the spectra for a given inflation- 
ary model (and were compared with the exact spectra of 
power-law inflation 0). 

The power spectra to next-to-leading WKB order, in 
the adiabatic expansion, are given by 



AD 

)S 



neim^yaHj (1 - ei) Wn, s 



(5) 



16 H 2 



'pi 



aH 



^ (l + gffx) 
(1 - ei) wii,t 



where mp\ is the Planck mass and all quantities are eval- 
uated in the super-horizon limit. A crucial part of our 
method is the evaluation of £jj and g^ ■ In particular, 
in this work we will estimate analytically next-to-leading 
WKB corrections to 0(ei) and show that the leading 
WKB order can also give predictions to (D(e 2 ). 

For the precise definition of gff and its evaluation we 
refer to Eq. (61b) and Section V of Ref. 8]. The quantity 



£77 was also defined in Eq. (26b) of the same reference, 
and can in general be written as 



£n(ffr,»7o) 



y/A 2 (n) -fc 2 77 2 — 
V 



(6) 



with 770 the time for which the integrand vanishes, r){ 
the super-horizon limit and A 2 (77) contains all the depen- 
dence on the horizon flow functions (77) . In a manner 
similar to repeated integration by parts, we can obtain 
general expressions valid for every A 2 {rf) which contain 
terms easy to evaluate explicitly and new integrals of suf- 
ficiently high order in the £4(77) so that they can be ne- 
glected. For example, on employing such a process once, 
we obtain the identity 



tniVuVo) = -\jA 2 {nf) - k 2 ril 



In 
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r ' a A'(rj) 



A(r)f) + y/A 2 { m )-k 2 rjf 



In 



'/f 



A(rj) - \jA 2 (n) - k 2 rf 



A(rf) + y/A 2 (n) ~ k 2 r/ 2 



dr?.(7) 



Since the integral in the right hand side contains A' , it is 
of (at least) one order higher in the ej than the remaining 
terms, which can be calculated explicitly. More details, 
omitted here for the sake of brevity, will be given in a 
forthcoming paper [Isj . 

Next-to-leading WKB order and first slow-roll order. The 
adiabatic corrections g^ in Eqs. (JSJ to leading order in 
the horizon flow functions are given by 



,AD 



37 19 



ff « s " 324 ~ 243 l ei+ 2 62 



1 



(8) 



f(l)T 
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324 ~ 243 



We can now write the expressions for the scalar and ten- 
sor spectra to next-to-leading WKB order (indicated by 
the subscript WKB*) and first slow-roll order (indicated 
by the superscript (1)) 



H 2 



C.WKB, " K ( xm 2^ — 



1 - 2 (D WKBt +l)ei- flwKB, e 2 - (2 ex + e a ) In ( — 



16 H 2 



A v 



/l, WKB* " nm 2 i -WKB, 



(9a) 



1-2 (D WKB , + 1) e x -2ei In ( — 



where A WKB * = 361/18 e 3 & 0.999 and -D W kb* = 777 — ln3 ~ —0.7302. The slow-roll approximation predicts, 
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for the corresponding quantities, A SR — 1 and D SR — 
C = 7 B + In 2 - 2 ps -0.7296 (where j E is the Euler- 
Mascheroni constant). Thus, the next-to-leading WKB 
order gives an error of about 0.1% for the estimate of the 
amplitude and one of about 0.08% on the coefficient C. 
This shows that analytical results obtained by the WKB 
method have reached the same accuracy as the standard 
slow-roll approximation. We also obtain the same slow- 
roll spectral indices and a-runnings, 

4!wkb, - 1 = -2 ei - e 2 , 4!Ikb. = ~ 2 £ i ( 9b ) 

4!Ikb, = ^WKB, = . ( 9c ) 
on using respectively Eqs. (|4b[) and l|4cfl. From Eq. (|4d|) 



the tensor-to-scalar ratio becomes 

R$L. = 16 e x (1 + D WKBt e 2 ) . (9d) 

Leading WKB order and second slow-roll order. Wc 
would now like to increase our accuracy in the slow-roll 
parameters to second order, while keeping the WKB ap- 
proximation to leading order (a further improved treat- 
ment to second order in the slow-roll parameters is in 
progress 0). On setting gff = in Eqs. J5J), we can 
write the expressions for the scalar and tensor perturba- 
tions to leading WKB order (indicated by the subscript 
WKB), and second slow-roll order (indicated with the 
superscript (2)) as 



,(2) 
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[-2 e x - e 2 + 2 (2 L» WKB + 1) ef + (2 D WKB - 1) e x e 2 + D WKB e 2 - D WKB e 2 e 3 ] In ( — 
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where A WKB = 18/e 3 ~ 0.896 and D WKB = | — In 3 
—0.7653. The spectral indices 114 bfl are then given by 



i (2) 

'S,WKB 



l (2) 
T,WKB 



-2ei-e 2 -2e 2 -(2D WKB +3) ex e 2 

" ^WKB 

e 2 e 3 (10b) 
-2ei-2e 2 -2 (D WKB + 1) e a e 2 



and their runnings (|4c|l by 



* (2) 

S.WKB 
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We note that the second order correction — 2 e\ to ns and 
nx agrees with the exact spectral index for power-law 



inflation 231. The tensor-to-scalar ratio 14dll becomes 



I? 



(2) 



= 16 ei 



1 + A, 



: e 2 



A, 
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Of course, Eqs. llOal) (|10c|) give the same results as ob- 
tained in Ref. Q , to first order in the horizon flow func- 
tions. As already found in Q, the leading WKB order 
gives an error of about 10% for the estimate of the ampli- 
tude and one of about 5% on C with respect to standard 
slow-roll results. The runnings as and «t are predicted 
to be of Q(e 2 ) vM - and in agreement with those obtained 
by the GFM [(J |2(| ■ It is notworthy that our power spec- 
tra and spectral indices have the same polynomial struc- 
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FIG. 1: Comparison, for the same inflationary model (with 
ei = £2 = 0.1, £3 = 0, k* = 0.05 Mpc -1 ), of angular power 
spectra for temperature anisotropies of the cosmic microwave 
background due to scalar perturbations with different accura- 
cies in the spectral slope ns, as in Eqs. lllObll and (llOct . The 
dotted line is the first-order slow-roll prediction, the dashed 
line has ns accurate to second order and no running, the dot- 
dashed has ns to first order and running as = —0.02, the 
solid line is the full second order prediction. The other pa- 
rameters used for this flat model are ficoivi = 0.26, fib = 0.04, 
Q.A = 0.7 and Ho = 72 km s -1 Mpc -1 . Let us note that the 
WMAP first year observational relative error on the first peak 
is 0.7% and on the second peak 1.2% plj. 



ture in the €i as those of the latter references. For the 



spectral indices ns and n T an analogous structure is also 
confirmed by the uniform approximation |jj 0, [n| [24| . 

Conclusions. We have shown that inflationary theoret- 
ical predictions have now reached expressions to sec- 
ond order in the slow-roll parameters confirmed by two 
completely different approximation schemes, such as the 
GFM and WKB. As can be seen from Fig. ^ f° r some 
inflationary models, second order slow-roll corrections to 
the power spectrum are necessary to perform a correct 
comparison between theoretical predictions and obser- 
vational data. The figure also shows that not only the 
runnings are important, but also the C(ef ) terms in the 
spectral indices. 

The different predictions of the two methods for the 
numerical coefficients in front of the O(e^) terms are at 
most of the order of 5% for the spectral indices and of 
10% in the amplitudes. For slow-roll parameters ~ 0.1, 
this leads to an accuracy of about 0.5% in the theoretical 
predictions for the tensor-to-scalar ratio llOdfl . A similar- 
accuracy (of 0.1%) has also been reached by the Boltz- 
mann codes p^. The predictions in the numerical coef- 
ficients and in the amplitude can be further improved by 
employing the next-to-leading WKB method [lg while 
including terms to second order in the slow-roll parame- 
ters. 
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